Bright and dark modes induced by graphene bubbles 
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Through a lattice dynamics analysis, it is revealed that the bubble plays a role of energy shield 
in the graphene, which helps to split the normal modes into two categories of distinct topological 
nature, namely the bright and dark modes. The topological invariants, Euler characteristic, of the 
bright and dark modes are 1 and 0, respectively. For bright modes, the energy is confined inside the 
bubble, so this type of modes are sensitive to the shape of the bubble; while opposite phenomenon 
is observed for the dark modes. The different behavior from these two types of normal modes is 
examined and verified in the process of phonon thermal transport. The bright and dark modes are 
expected to be distinguished in experiment with existing scanning force microscope techniques, and 
they should play significant roles in many other physical processes. 
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The strain engineering provides an exciting approach 
to improve or manipulate various properties of graphene- 
based materialsii By generating strain, the experimen- 
taUsts have been able to study the zero-field quantum 
hall effect in graphene and realize energy gap with the 
help of the stained supperlatticeS)^ and observe the strain 
induced red-shift of the G modei^i^ The strain effect 
has also been theoretically investigated for different as- 
pects of graphene, such as the strain induced electro- 
chemical potential and charges^^ the electronic band 
gap for single layer graphene^ and bilayer graphene^, 
the pseudomagnetic fields, ^"^^ the optical properties, 
and the hydrogen storage capabilities of metal-decorated 
grapheneJ^ As a specific morphology of strain, a bub- 
ble in graphene was observed by Stolyarova et.al in 
2009, which is formed due to the gas between graphene 
and its substrate.— In 2010, Levy et.al found a very 
strong pseudo-magnetic field inside the bubble by us- 
ing the scanning tunneling microscopy to do the spectro- 
scopic measurements Quite recently, the Manchester's 
graphene research group show that it is possible to con- 
trol the shape of the bubble in graphene by an external 
electric field i'' This technique is proposed to produce op- 
tical lenses with variable focal length. 

In this paper, we investigate the effect of the bubble 
on the normal modes in graphene. We distinguish two 
opposite categories of normal modes, namely bright and 
dark modes, which are topologically different. The bub- 
ble is full of energy for the bright modes, while there 
is no energy inside the bubble in the dark mode. Due 
to their topological difference, the bright/dark mode is 
sensitive/insensitive to the shape of the bubble in the 
graphene. We examine the different behavior of bright 
and dark modes in the phonon thermal transport pro- 
cess. It is expected that the distinction between these 
two types of normal modes should be able to be mea- 
sured in the experiment. 

Figure [1] displays the configuration of a bubble within 
the graphene. The size of the graphene is 59.6 Ax44.3 A. 




FIG. 1: Two different views of a bubble in the graphene. The 
size parameters of the bubble are {h, A) = (7, 8) A, where h is 
the height of the bubble and A reflects the lateral size of the 
bubble. 



The xy-plane lies in the graphene plane and z-direction 
is perpendicular to the graphene plane. The bubble has 
a Gaussian shape. The z-coordinate of atoms on the 
surface of the bubble is given by z{x, y) = /i*exp{— [(a; — 
xqY + {y~yaf]/{2X^)}, where (a;o,yo) are the x- and y- 
coordinates of the center of the bubble. The bubble thus 
is characterized by the height h and a lateral dimension 
A. (/i. A) = (7, 8) A for the specific bubble shown in the 
figure. 

We calculate the normal modes by diagonalizing dy- 
namical matrix, which is obtained from the spring force 
constant model with longitudinal parameter ki and trans- 
verse parameter k±_ upto the fourth- nearest neighbors)^ 
The valence force field model^^ has also been tried and 
the calculation results are qualitatively the same. Calcu- 
lations shown below are based on the spring force con- 
stant model. For each normal mode (w,^), the energy 
carried by atom j is proportional to ^jj'^iS.'jx + + S.'jz)^ 
where with a = x,y, z arc the three components cor- 
responding to atom j in the polarization vector ^. We 
analysis the energy spatial distribution for each normal 
mode. Fig. [5] shows the bubble energy in logscale for all 
normal modes in (a) pure graphene, (b) graphene with 
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FIG. 2: The energy within bubble regions for all normal modes in (a), pure graphene; (b). graphene with bubble (7, 8) A; and 
(c). graphene with bubble (15, 8) A. The y axis is plotted in logscale. 
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FIG. 3: The energy spatial distribution of bright modes in the top four panels, and of dark modes in the bottom four panels 
for graphene with bubble (7, 8) A. Value in each panel is the frequency of that normal mode. 



bubble (7, 8) A, and (c) graphene with bubble (15, 8) A. 
The bubble energy of a normal mode is defined to be the 
energy of this mode confined within a circular region of 
radius 5 A around the bubble center. The bubble energy 
has been normalized in such a way that the total energy 
of the normal mode is unitary. Panel (a) shows that the 
bubble energy of the pure graphene is about 3% for all 
normal modes. It indicates that the energy is uniformly 
distributed on the entire graphene, considering the area 
of the bubble region is also around 3% of whole area of 
the graphene. Panel (b) exhibits strong influence from 
the bubble on the energy spatial distribution of normal 
modes in graphene. On the one hand, the bubble en- 
ergy of many normal modes is distinctly enhanced, and 



can reach as high as 80% for some normal modes. On 
the other hand, the bubble energy of some other normal 
modes are strongly suppressed to be two or three orders 
smaller than the pure graphene's value of 3%. These 
observations manifest that the bubble acts as an energy 
shield. It can thus split the normal modes in graphene 
into two opposite categories. One is to confine the energy 
inside the bubble region, while the other tries to repel en- 
ergy for the bubble region. Panel (c) further reveals that 
the bubble energy is more seriously enhanced/suppressed 
by taller bubbles; thus they show a better performance 
as an energy shield. 

The energy spatial distributions are shown in Fig. |3] for 
eight normal modes in graphene with bubble (7, 8) A. 
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FIG. 5: The circular area (left) with the Euler characteristic 
of 1 and the cirque area (right) with the Euler characteristic 
of 0. 



FIG. 4: (Color online) The histogram of bright and dark 
modes in graphene with bubble (7, 8) A. The total numbers 
of the bright and dark modes in this system are 101 and 112, 
respectively. 



theorem, after applying to these two special geometries: 
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The four upper normal modes belong to the category 
with large bubble energy. The energy is concentrated 
so strongly in the bubble region that the bubble looks 
bright. For this reason, we name this type of normal 
modes as 'bright modes'. The opposite behavior is ob- 
served in the four lower panels. The bubble energy is 
almost zero, which means that the bubble region is an 
energy vacancy in these modes. The bubble looks dark 
and similarly this type of normal modes are named 'dark 
modes'. The current scanning force microscope tech- 
nique is able to image the vibrational morphology of 
normal modes, such as the edge modes in the graphene 
nanoribbon.— Hence, the bright and dark modes we have 
found here are readily to be observed in the experiment. 
The distinctness between these two type of normal modes 
is also quite possible to be verified experimentally. Fig. |4] 
locates the position of bright and dark modes in the fre- 
quency space for graphene with bubble (7, 8) A. It shows 
that the bright mode is more uniformly distributed in 
whole frequency range while dark modes mainly locates 
in the high frequency region. This figure indicates that 
bright mode will play a role in whole temperature range. 
However, the dark mode will become important at high 
temperatures, where high frequency normal modes can 
be sufficiently excited. 

From the topological point of view, we can clearly dis- 
tinguish the bright and dark modes. The topological na- 
ture of the bright and dark modes is equivalent to the 
circular and cirque respectively as shown in Fig. [5] The 
gray filled area corresponds to the high density of spatial 
energy for normal modes. The Euler characteristic is the 
topological invariant. It coincides with the Euler num- 
ber for the two geometries in Fig. [5] The Euler charac- 
teristic (x) can be calculated through the Gauss-Bonnet 



where Kg is the curvature and ds is the line element 
along the edge. Using this formula, we get the Euler 
characteristic for the circular and the cirque are 1 and 
0, respectively. Hence, the Euler characteristic for the 
bright mode is 1 and that of the dark mode is 0. 

From the above, we have learn the energy spatial dis- 
tribution for the bright and dark modes, from which we 
demonstrate that these two types of normal modes are of 
different topological nature with different Euler charac- 
teristic. As we know, the topological invariant is of key 
importance in many physical processes. It is quite possi- 
ble that the different topological invariant of the bright 
and dark modes will lead to observable distinctness in 
those physical processes where the normal mode is an im- 
portant exciton. The phonon thermal transport is one of 
such physical process, where normal mode is the energy 
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FIG. 6: (Color online) The transmission for bright and dark 
modes versus the shape of the bubble. The two straight lines 
are the fitting curves for two data sets. 
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FIG. 7: The thermal conductance of graphene with bubbles 
of varying shape at room temperature. The straight line is 
guide to the eye. 



carrier. The thermal transport in graphene is currently 
an active field, where very high thermal conductivity was 
found 1^"— For a recent review of this subject, we refer 
to Ref. [H. We apply NEGF approach^ to calculate the 
transmission function, T[uj], of normal modes. 



T[lo] = Tr (GTiG'^r/ 



(2) 



where G° = (C)^ is the advanced Green's function 
and r^/zj are the two self energy. The phonon ther- 
mal conductance is calculated by the Landauer formula. 
Fig. ini shows the transmission for a bright mode at 
ui = 792cm~^ and a dark mode at = 1622cm~^ in 
graphene with bubbles of different shapes. Obviously, 



the transmission for the bright mode is very sensitive to 
the geometrical ratio h/X, and decreases linearly with 
the increase of h/X. It indicates that this mode becomes 
more bright for a bubble of lager h/X, where the energy 
is more concentrated in the bubble region. On the con- 
trary, we find that the dark mode is insensitive to the 
ratio of h/X. This is because the energy of the dark 
mode is distributed throughout the whole graphene sam- 
ple except the bubble region. Thus it is insensitive to the 
local environment of the bubble. The room-temperature 
thermal conductance result in Fig. [7] confirms that the 
bright mode plays important role and leads to a linearly 
decreasing of the thermal conductance with the increase 
of h/X. 

In conclusion, we have studied the normal modes in 
graphene with bubbles and found two topologically dif- 
ferent types of normal modes, namely the bright and dark 
modes. The bright modes confine energy inside the bub- 
ble region and has an Euler characteristic of 1; while 
the dark mode repel energy for bubble region and has 
an Euler characteristic of 0. As a result of their differ- 
ent topological nature, the bright mode is sensitive to 
the shape of the bubble; yet the dark mode is insensitive 
to the shape of the bubble. The distinctness between 
bright and dark modes is expected to be observed in the 
experiment. We have also examined the difference be- 
tween these two types of normal modes in the phonon 
thermal transport process, and find that the transmis- 
sion of bright/dark mode is sensitive/insensitive to the 
geometrical ratio h/X. 
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